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' As a typical representation of complex networks studied relatively thoroughly, financial market 

\ presents some special details, such as its nonconservation and opinions spreading. In this model, 

. agents congregate to form some clusters, which may grow or collapse with the evolution of the 

system. To mimic an open market, we allow some ones participate in or exit the market suggesting 
that the number of the agents would fluctuate. Simulation results show that the large events are 
frequent in the fluctuations of the stock price generated by the artificial stock market when compared 
pH ■ with a normal process and the price return distribution is a levy distribution in the central part 

followed by an approximately exponential truncation. 
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Ch ■ I. INTRODUCTION 
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As the research of the complex systems is getting deeper and deeper, to find the universal rules and principles of 
these systems and to answer the origination of "complexity" become more and more attractive 0, 0, 0, SIE 0| ■ Thus, 
the vision of physicists is no longer restricted in the traditional areas but concentrates on the more comprehensive 
domains, leading to the birth of many burgeoning disciplines through the interaction and amalgamation of physics 
' and other fields such as biology, finance, sociology, and so on. Over the last decade, there has been significant 
' interest in applying physical methods in social and economical science . In particular, the study of financial market 
04 prices has been found to exhibit some universal properties similar to those observed in physical systems with a large 
number of interacting units, and some models m iH III (as we know, the first stock market simulation was 
performed by economist Stigler in 1964 j^J) have been introduced to the financial and more recently to the physical 
communities which attempt to capture the complex behaviors of stock market prices and market agents Il4l . These 
models, including behavior- mind model 15^ J^], dynamic-games model J/^], multi-agent model 18. 19i l20ll2lLl23 .l2 ^ 
and so on, are based on the statistical properties of price fluctuations which should be recovered by more suitable 
microscopic models |2^: (l)sharp peak and fat-tail distribution for the price changes (the return histogram) [26l| : 
(2)the distribution of returns decays with power law in the tails |27j . with exponent near 3; (3)price fluctuations are 
^ ] not invariant against time reversal, i.e. they show a forward-backward asymmetry p^ . 

O ■ Among the more sophisticated approaches are dynamic multi-agent models based on the interactions of two dis- 
^ tinct agent populations ("noisy" and "fundamental" traders), which could simulate the price forming processes and 
^ reproduce some of the stylized observations of real markets, but fail to account for the origin of the universal char- 
acteristics. An alternative approach, the herd behavior .30] explored in this paper, may be capable to induce the 
power-law asymptotic behavior in the tail of price return distribution with an exponent well fitting the truncated Levy 
regime [2^ as found in real data 0, |23| • This approach has been formalized by Cont and Bouchaud ^| as a static 
percolation model. Subsequently, this percolation model has been bettered by introducing a feedback mechanism 
between the price return Z and trader activity a: a — > a -I- aZ , where a is the factor denoting the sensitivity to price 
fluctuations. Then the volatility clustering can be reproduced, and all of the statistical properties of fluctuations for 
prices mentioned above could be observed. 

In the Cont-Bonchaud model, random percolation clusters are used as groups of traders. In the simple version, at 
each iteration each cluster buys with probability a, sells with probability a or sleeps with probability 1 — 2a. When 
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FIG. 1: (a)To illustrate, this is a small-scale matrix comparing to our model which can help us to explain what is a cluster 
called in this paper. A cluster, here we call it M, is defined as an aggregation of nodes that share the same information and 
hold the same opinion. Moreover, in the topology, every two nodes in M can reach each other. Obviously, there are two 
different clusters which are represented by different colors, yellow and blue in figure 1(a). In the process of growth, each cluster 
has distinct probability to absorb new members who could occupy the frontier and empty nodes with different probabilities 
determined by the cluster which they would participate in. And the sites which the new agents could occupy are denoted by the 
hollow nodes. But the colors of the hollow nodes are distinguished which implies that there must be some difference between 
the positions marked by different colors. In fact, the positions denoted by red hollow nodes mean that if a new member comes 
into the market and takes up one of these positions, he would be puzzled for there are two different attitudes. In other words, 
he comes across two different information resources which would lead to completely different consequences. And he operates 
just like a information bridge which would stimulate the communication of clusters holding different attitudes and even lead 
to a serious confliction that could promote some certain clusters to merge others and expand to a giant one such as the cluster 
shown in figure 1(b). The opposite process called "collapse of the cluster" in this model, which means that the nodes belonging 
to a same cluster are removed from the lattice, reflects another phenomenon of the real market that there are always some 
traders quitting the market due to different reasons which may lead to a serious influence on the price. 

the activity a is small, there are only few clusters trading at a time most of the times. Therefore, the distribution 
Pb{Z) of relative price fluctuations or "returns" Z scales as the well-known [s^, H^i 01 cluster size distribution of 
percolation theory. But when we increase the value of a , more and more clusters would make contribution to the 
relative price fluctuations and the central limit theorem suggests that the distribution Pb{Z) convert from power- 
law tails to a more Gaussian shape for large systems. Price changes in the logarithm (return) are proportional to 
the difference between the supply and demand. On average, price rises or falls with equal probability and without 
correlations between consecutive steps. An assumption made in this model which should not be ignored is that the 
probability a (activity) is set to be the same for all groups and remains constant through out the whole process, which 
may be a good strategy for simplifying a physical model but may be not a good regulation for establishing a model 
which we expect to reflect the various phenomena found in the real stock market as genuine as possible so it could be 
more helpful for us to capture the complex properties of the real world. Although in the successional studies of the 
percolation model different mechanisms are used to establish a self-organized model where the investor groups with 
various trading activities and sizes are formed automatically, few of these models considered the fluctuations of the 
traders in the market, fact that there are always agents who take part in or exit from the market due to various reasons 
which might have a serious influence on the price. Here we introduce a self-organized model where the activities and 
sizes of different investor groups are driven by the confliction and harmonization of the strategies adopted by different 
groups. The simulation results which well agree with the observations of real markets are also shown in the third 
section. 



II. THE NEW MODEL 

Considering an open market which absorbs and removes traders with probabilities respectively in the process 
of trading, we draw the growth and collapse of clusters into our model based on Cont-Bouchaud model's random 
percolation. After every N iterations of trading, each cluster which is defined in figure 1 grows around itself (figure 
1(a)) with probability P^, collapses and annihilates with probability Pm or sleeps with probability 1 — Pd — Pn- Once 
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a cluster collapses, some new clusters will come into the market at the positions where the collapse occurred (figure 
1(b)), with a fixed probability P^. We build our model as follows based on the thoughts above: 

(1) Initiation: a, L x L lattice is occupied randomly with probability Pin, and each cluster is randomly given a state: 
buying, selling, or sleeping, which are represented by 1, -1 and respectively. 

(2) Trading: at time t, each agent in the market sells or buys a unit of stock. Then we calculate the difference 
between the supply and demand 



(1) 



where m is the total number of agents who are presented in the market at t and Si represents the state of the i-th 
agent. The evolution of price follows the rule: 

P(i + 1) = P(t)e'^/^ (2) 

where X — Np + Nn- Here Np and A''„ denote the number of buyers and sellers, respectively. Figure 2 shows the price 
time series, which is rather similar to that of real stock market. And the active probability 2a with which the agents 
choose buying or selling rather than sleeping evolves following the Equ.(3): 

a{t) = a{t -l) + lr (3) 

where I represents the sensitivity of activity a{t) to the difference between the demand and supply. And then, each 
cluster buys, sells, or sleeps with probabilities 2a{t)pi,, 2a{t)ps, 1 — 2a{t) respectively, where 

r /i + lyir, r<0 , . 

+ V2r, r > ^ > 

and 

Ps = l- Pb (5) 

The first term on the right hand of both of Equ.(4) denotes the probability with which the active one would buy 
rather than sell without considering the feedback of the price fluctuations. The difference between the case r < and 
r > is the coefficient of the last term on the right hand and takes into account that agents are risk adverse and 
thus more impressed by a downturn than by an upturn of the market so that the parameters vi and 1^2 denote the 
sensitivity of the agent's mentality to the price fluctuations. (The price would fall if r < and would rise if r > 
according to Equ.(2) 35, 36, 37, 38j.3£, j|£|). So the value we adopt of i^i is smaller than that of 1^2 |4ll |. 

(3) Growth: after each N iterations, there are three types of evolution with different probabilities respectively 
depicted in the following segments: growing, collapsing and sleeping. 

The first situation is that new traders come into the market occupying the empty sites around the old clusters, for 
example, cluster G, just as the sites marked by hollow nodes proposed in figure 1(a) with the probability 

P!{t + l)=P!{t) + k{NT-c^{t)) (6) 

where fc is a kinetic coefficient and Nt is a threshold parameter and 

nig 

in which rUg represents the scale of the cluster G, in other words, the number of nodes which belong to G. s| denotes 
the state of the j-th node belonging to G. The probability Pj'(t + 1) is obviously limited to the range [0, 1] so that 
we have to impose Pj'(t + 1) = and = 1 if the recurrence relationship Equ.(6) gives values for Pj'(< + 1) < or > 1. 
If a few of different clusters whose states are different encounter, one ( noted by V ) will defeat others ( the total 
number is n, including V ) with the probability 

n 

P.{t)^\c^t)\/Y,\c\t)\ (8) 



And the evolution due to Equ.(8) would lead to the consequence that the defeated clusters would accept the opinion 
and adopt the same strategy of the winner. In other words, they are annexed by the winner. By contrary, when the 
states of the encountered clusters are all the same, they would combine and operate as a whole. 
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FIG. 2: Time series of the typical evolution of the stock price, where Ph=0-01. One can see that the trend and fluctuations of 
the stock price generated by our model are rather similar to that of real stock market. 
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FIG. 3: The returns of simulated price fluctuations for St=l. It can be seen that large events are frequent in the fluctuations 
of stock market when compared with a normal process. 



The second, each old cluster such as Q collapses with the probability 



(9) 



which indicates that the probability with which a cluster collapses would increase with its growth. When a cluster 
takes up all of the sites of the lattice, it would surely collapse. Once the old cluster collapses, the members of the new 
clusters whose states arc not necessarily the same as the old one are produced with fixed probability Ph and take up 
the sites where the old cluster has existed. 

The final circumstance, clusters sleep with the probability 1 — Pd{t) — Pn{t)- 

(4) Repeat step (2) and (3) for enough times. 
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FIG. 4: The probability distributions of price returns with (5f=l,2,4,8,16,32,64 respectively. In this figure, the central part of 
the distribution of returns appears to be well fitted by a levy distribution. 



III. SIMULATION RESULTS 



The typical parameter space we adopt in our simulation is as follows: I — 0.0001, k ~ 0.0001, Nt — 50, Ph — 0.01, 
a(0) = 0.35, Pd{0) = 1.0, N = 100, /i = 0.59, i^i = 0.00005 and 1^2 = 0.0001. About 100 traders (in other words, 
^'in=0.01) are distributed on a square lattice randomly and the initial stock price is 1.0. The simulation results are 
very sensitive to some of the parameters such as I, k, fi, vi and V2. When the values of them are little larger, the 
price fluctuations would be very exquisite and when they are little smaller, the price trend would be very meek. But 
the simulation results are not very sensitive to other parameters for other values of them could lead to the results 
which are in good agreement with the real data, too. But there are some amazing results we should point out: (1) /i 
is not 0.5 but 0.59, little larger than 0.5, which is very close to the threshold value: 0.593 in the percolation theory 
which may be a support to the point that the real markets should operate close to the critical point where profitable 
trade opportunities are barely detectable The process by which the market self-organizes close to the critical 

point is more likely to be of evolutionary nature and hence to take place on longer time scales '43|. And this result 
suggests that the choice whether to buy or to sell is not completely random as the traditional point stands which 
implies that should be 0.5. (2) From the value space we could see that V2 = 2vi, which means that the affliction 
which is brought by losing one unit of wealth would be twice as much as the satisfaction which is brought by gaining 
the same amount of wealth according to Kahneman's Prospect Theory 40] . 

To compare the statistical properties of the price generated by our model and that of real stock markets further, 
we study the returns of price which is defined as Equ.(lO): 

Z {t) ^ log P{t + 5t)~ log P{t) (10) 



Mandelbrot has proposed that the distribution of returns is consistent with a Levy stable distribution 45]. In 1995, 
Mantegna and Stanley analyzed a large set of data of the S&P500 index. It has been reported that the central part 
of the distribution of S&P500 returns appears to be well fitted by a Levy distribution, but the asymptotic behavior of 
the distribution shows faster decay than that predicted by a Levy distribution 26, 46]. The similar characteristic of 
the distribution of returns is also found in Hang Seng index • Figure 4 shows the probability density of normalized 
returns, which display a clear Levy distribution for 5t=l, 2, 4, 8, 16, 32, 64. 

Because larger 5t implies less data points, it is difficult to determine the parameters characterizing the distributions 
only by investigating the spreads. Hence, we studied the peak values at the center of the distributions, i.e., the 
probability of zero return Pb{Z = 0) as the function of 5t. With this choice, we can investigate the point of each 
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FIG. 5: The central peak value as a function of St. The slope of the fitted line is —0.61 ± 0.01 which is very close to the real 
value 0.62 found in Hang Seng index [i^. 




FIG. 6: Re-scaled plot of the probability distributions shown in figure 4. Data collapse is evident after using re-scaled variables 
with Q — 1.61. The abscissa is the re-scaled returns, the ordinate is the logarithm of re-scaled probability. 



probability distribution which is least affected by noise. Figure 5 shows the central peak value versus St in a log-log 
plot. It can be seen that all the data is well fitted by a straight line with the slope —0.61 ± 0.01 which is very close to 
the real value 0.62 found in Hang Seng index This observation agrees with theoretical model leading to a Levy 
distribution. 

If we assume that the central part of the distribution of returns can be described by a Levy stable symmetrical 
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FIG. 7: The accumulate probability distributions Pb{g > Z) of 1-minute returns generated by our model. For data in the 
region 10 < Z < 200, regression fits yield a = 2.93 (positive tail) and a = 2.78 (negative tail). 

distribution with an index a and parameter 7, 

/•oo 

P;^{Z,St) = (I/tt) e-T"'*l9l° cos(gZ)d(7 (11) 
"'0 

where e""*""^*''' is the characteristic function of a Levy symmetrical stable process, the probability of zero return is 
given by 

n(0) = Pr(0, St) = r(l/a)/[™(7<50'/"] (12) 

where F is the Gamma function. Using the value —0.61 ± 0.01 for the slope of the fitted line to the data (figure 5), 
we obtain the index a = 1.61 ± 0.02. 

To check whether the Levy scaling can be extended to the entire probability distribution of returns generated by 
our model, we notice that under the transformation: 

Z, = Z/[(<5i)(i/")] (13) 

and 

Pis{Zs) = {6tY'''Pb{Z,5t) = {5tY'''P^[{5tY'''Zs,5t] (14) 

the distributions for different time scales 5t will collapse onto one curve. Figure 6 shows the re-scaled distributions 
for the same data in figure 4 in the scaled variables, i.e., PbsiZg) versus Zg. Data collapse is evident, except for some 
data points in the tails for large St. The closer to the central point Zg = 0, the stronger is the extent of data collapse. 
These observations imply that the Levy distribution is a better description of the dynamics of the random process 
underlying the variation of returns in the central part of the probability distribution Pb{Z) over St spanning at least 
two orders of magnitude. 

In order to determine if an exponential truncated Levy distribution can be used to describe the stochastic process 
and to investigate the kind of asymptotic behavior outside the Levy stable region, we study the accumulate distribution 
Pb{g > Z) of the fluctuations of financial data. 

For a stable symmetric Levy distribution (0 < a < 2), the two tails show a power-law asymptotic behavior 



(15) 
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and hence the second moment diverges. This leads to an asymptotic power-law for the accumulate distribution for 
both the positive and negative tails 45] in the form 

Pi,{g >Z)^ (16) 

Figure 7 shows the accumulate probability distribution of returns Ph{g > Z) for 6t = 1 min. for the data generated 
by our model. For data in the region 10 < Z < 200, regression fits yield a = 2.93 (positive tail) and a = 2.78 
(negative tail). These results appear to be outside the Levy stable range of < a < 2 but they fit well the result 
produced from the real data which is near 3. 



IV. CONCLUSION 

Compared with the Cond-Bouchaud percolation model, our model presents a nonconservation market. With the 
evolution of the model's topology, there are new traders coming in and old ones leaving, which depicts the real stock 
markets more approximately. What is more, the process of the amalgamation and expansion, and the breakdown of 
the clusters in our model well consists with the phenomena (so called herd behaviors) in the real stock market that 
more and more people would take the same performance when they found more and more people around them take 
the same action. Some other simulations show that the information entropy, when we consider the clusters as various 
information resources and the process of merging and collapsing as the spreading and dying out of the information, 
has some relationships with the point which we define as the break-point of price. Moreover, there are also some 
amazing facts in the results of our simulations such as why the first term on the right hand of Equ.(4) approximate 
to the threshold value in the percolation theory. Since the main goal of this article is to establish and describe the 
model itself, we would not give detailed simulation results and analysis, which will be given elsewhere soon. 
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